Astrophysical factor for the neutron generator 13 C(a,n) ie O reaction in the AGB stars. 
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The reaction 13 C(cv,n) is considered to be the main source of neutrons for the s-process in AGB 
stars. At low energies the cross section is dominated by the l/2 + 6.356 MeV sub-threshold resonance 
in 17 whose contribution is determined with a very large uncertainty of ~1000% at stellar temper- 
atures. In this work we performed the most precise determination of the low-energy astrophysical 
S factor using the indirect asymptotic normalization (ANC) technique. The a-particle ANC for the 
sub-threshold state has been measured using the sub-Coulomb a-transfer reaction ( 6 Li,d). Using 
the determined ANC we calculated S(0), which turns out to be an order of magnitude smaller than 
in the NACRE compilation. 

PACS numbers: 25.70.Hi, 25.55.Hp, 26.20.+I, 27.20.+n 
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About half of all elements heavier than Iron are pro- 
duced in a stellar environment through the s-process, 
which involves a series of subsequent neutron captures 
and /3-decays. The reaction 13 C(a,n) 16 is considered to 
be the main source of neutrons for the s-process at low 
temperatures in low mass stars at the Asymptotic Gi- 
ant Branch (AGB) 0. This is because it is exothermic 
and can be activated at low temperatures. Two factors 
determine the efficiency of this reaction: the abundance 
of 13 C, and the rate of the 13 C(a,n) reaction. Accurate 
knowledge of the 13 C(a,n) 16 reaction rates at relevant 
temperatures (0.8 - 1. Ox 10 s K) would eliminate an es- 
sential uncertainty regarding the overall neutron balance 
and will allow for better tests of modern Stellar models 
with respect to 13 C production in AGB stars (see Q and 
references therein). 

The rate of the 13 C(a,n) reaction at temperatures of 
~10 8 K is uncertain by ^300% 0] due to the prohibitively 
small reaction cross section at energies below 300 keV. A 
directly measured 13 C(a,n) cross section is only avail- 
able at energies above 279 keV (see and references 
therein). Below this energy the cross section has to be 
extrapolated. It was shown 3, 4] that this extrapolation 
can be strongly affected by the l/2 + sub-threshold res- 
onance in 17 at 6.356 MeV excitation energy, which is 
just 3 keV below the a threshold. It was assumed in the 
recent NACRE compilation Q that this resonance has 
a well developed a cluster structure. This assumption 
leads to a strong enhancement of the cross section at low 
energies ||. Recently Kubono et al. Q determined the 
contribution of the sub-threshold state at 6.356 MeV in 
17 to the astrophysical factor for the 13 C(a,n) reaction 
at low energies by measuring the a-particle spectroscopic 
factor of this state by the a-transfer reaction 13 C( 6 Li,d) 
at 60 MeV. The extracted spectroscopic factor was found 



to be very small, S a ~ 0.011 5], making the influence of 
this sub-threshold state on the astrophysical factor negli- 
gible. However, it was shown in [|| that the same experi- 
mental data was compatible with a large S a factor for the 
sub-threshold state in question. It is the main goal of this 
work to resolve this difference and to develop a technique 
which determines the contribution of sub-threshold reso- 
nances to the (a,n) reaction cross sections using a model- 
independent approach. Until now the ANC method has 
been applied to determine the astrophysical factors for 
radiative capture processes 0, H> El • Here we present the 
first case of application of the ANC method to determine 
the astrophysical factor for the 13 C(a,n) 16 reaction. 

The amplitude of the reaction x+A — > b+B proceeding 
through the sub-threshold resonance F is given in the R- 
matrix approach by [Tol ] 
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Coulomb modified Whittaker function, ro the channel ra- 
dius, C£ A — C^ A /T(l + I + 77) stands for the Coulomb 
modified ANC for the virtual decay (synthesis) F «-> 
x + A, rj and I are the Coulomb parameter and rel- 
ative orbital angular momentum of the sub-threshold 
bound state (xA), and Tf(EbB, 7"o) is the resonance 
width for the decay to the final channel b+B. We assume 
that the total width of the resonance F is equal to Tf, 
Eij = kfj/(2 Hij) is the relative kinetic energy of particles 
i and j, K = \J2 [i xA £ and Ef is the binding energy for the 
virtual decay F — > x + A. In this case Tf = F„ = 124±12 
keV is a known neutron partial width. Thus, the 
ANC is the only missing quantity needed to calculate 
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Figure 1: Spectrum of deuterons from the 6 Li( 13 C,d) reac- 
tion, at 6° with the 8.5 MeV 13 C beam. The inset shows the 
level scheme of 17 The solid line is a Gaussian fit. 



the cross section for the 13 C(a,n) 16 reaction proceed- 
ing through the sub-threshold state. Due to the periph- 
crality of the sub-Coulomb transfer reactions the over- 
all normalization of the a-transfer reaction cross section 
is determined by the product of the squares of the ini- 
tial and final ANCs rather than the spectroscopic fac- 
tors. The initial ANC for the a + d — > 6 Li is known, 
(Cl L j) 2 = 5.3 ±0.5 fm- 1 [12. Hence, by normalizing the 
DWBA cross section to the experimental one we can de- 
termine the ANC for a + 13 C -> 17 (6.356 MeV, 1/2+). 

In this Letter, we report the application of the ANC 
technique to determine the astrophysical S factor of the 
13 C(a,n) 16 reaction at astrophysically relevant ener- 
gies by measuring the ANC for the virtual synthesis 
a + 13 C -> 17 (6.356 MeV, 1/2+) using the a-transfer 
reaction 6 Li( 13 C,d), performed at two sub-Coulomb en- 
ergies, 8.0 and 8.5 MeV, of 13 C at the Florida State Uni- 
versity Tandem-LINAC facility. The choice of inverse 
kinematics, 13 C beam and 6 Li target, allowed measure- 
ments to be made at very low energies in the cm., and to 
avoid background associated with the admixture of 12 C 
in the 13 C target. Angular distributions of the deuterons 
from transfer reactions at sub-Coulomb energies in in- 
verse kinematics peak at forward angles. Four Si AE — E 
telescopes were positioned at forward angles to identify 
deuterons. Thicknesses of the AE detectors were in the 
range from 15 to 25 fim. 50 fig/cm 2 Li targets (enriched 
to 98 % of 6 Li) were prepared and transported into a 
scattering chamber under vacuum to prevent oxidation. 
The telescope at the smallest angle (6° in Lab. frame) 
was shielded from the Rutherford scattering of 13 C on Li 
target with a 5 /im Havar foil. A spectrum of deuterons 
at 6° (which corresponds to 169° in the cm. for the l/2 + 




cm. 

Figure 2: Angular distribution of the 13 C( 6 Li, d) 17 0(l/2+) 
reaction. Data taken at 8.5 MeV of 13 C are shown as dia- 
monds, and 8.0 MeV data as boxes. Dashed and solid lines 
are DWBA calculations at 13 C energies 8.31 and 7.81 MeV 
(see text). 



6.356 MeV state) at a beam energy 8.5 MeV is shown in 
Fig. ^ The typical experimental resolution in the cm. 
system (mainly defined by the 380 keV energy loss of the 
13 C beam in the 6 Li target) was about 250 keV (FWHM). 
The angular distributions of 6 Li( 13 C,d) 17 0(l/2+, 6.356 
MeV) are shown in Fig. [21 Their shape is typical for sub- 
Coulomb transfer reactions. Absolute normalization of 
the cross section was performed by measuring the elastic 
scattering of 6.868 MeV protons on the 6 Li target. The 
cross section of this reaction at 95° is known with 3% 
accuracy 0. Each telescope was sequentially placed 
at 95° and the product of the target thickness times 
the telescope solid angle (t x Af2) was determined for 
each telescope. The code FRESCO (version FRXY.3h) 
[j"H was used to calculate the angular distribution of the 
13 C( 6 Li,d) reaction in the Distorted Wave Born Approx- 
imation (DWBA) approach. The DWBA calculations 
were performed for beam energies at the center of the tar- 
get, 8.31 MeV for 8.50 MeV measurements and 7.81 MeV 
for 8.0 MeV energy of the beam. It was found that nor- 
malization factor is the same for both energies, indicating 
that Compound Nucleus mechanism plays only minor (if 
any) role. The extracted ANC, unlike the spectroscopic 
factor, does not depend on the number of nodes of the 
a — 13 C!/ 2 + bound state wave function, or the geometri- 
cal parameters of the Woods-Saxon potential. Parame- 
ters of optical model potentials of the usual Woods-Saxon 
form, used in the DWBA calculations, are given in Table 
HI The LCI potential was used for the 6 Li+ 13 C channel. 
It reproduces experimental data on elastic scattering of 
6 Li by 13 C at energies ranging from 3 to 23 MeV in cm. 
Experimental data on elastic scattering of deuterons on 
17 at low energies are not available. Thus, several po- 
tentials for the d + 17 channel were used \mUM An- 
gular distributions shown in Fig. [2 were calculated with 
the DPI potential, however it was verified that other po- 
tentials [lalla produce essentially identical results, with 
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Figure 3: S-factor of the 13 C(a,n) reaction. The experimental 
data from direct measurements, corrected for electron screen- 
ing, are from The contribution of l/2 + state is shown 
as the dashed curve. The dash-dotted curves represent 26% 
uncertainty band. 
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Figure 4: Rate of the 13 C(a,n) reaction. The dash-dotted 
curve is from solid curve is the rate obtained in this work. 
Yellow band in inset represents uncertainties from |3j, new 
uncertainties are shown in red. 



variations in normalization factor of less than 7%. This 
demonstrates that the transfer reaction cross section at 
sub-Coulomb energy only weakly depends on the param- 
eters of the optical potentials, which was the main point 
of making this measurement. In fact, calculation with 
no nuclear part of the optical potentials changes the ab- 
solute value of the cross section at large angles by only 
^40%. Further investigation of the cross section sen- 
sitivity to the parameters of the optical potentials was 
performed and was found to be less than 20% if param- 
eters kept within reasonable limits. We found no sen- 
sitivity of the extracted ANC to the parameters of the 
core-core DC1 interaction potential in the full DWBA 
transition operator. Our determined Coulomb- modified 
ANC squared for 13 C + a -> 17 0(l/2+,6.356 MeV) is 

(^r°c /2+) ) 2 = 0.89 ± 0.23 far 1 . The contribution of 
the l/2 + state to the astrophysical S factor calculated 
using Eq. (JIJ is shown as a dashed curve in Figure It 
was verified that this result is insensitive to variations of 
the channel radius. 

Five sources of uncertainty associated with the S(0) 
factor of the 6.356 MeV l/2 + state (same as for the 
ANC), can be identified: 7% statistical uncertainty, 7% 
combined systematical uncertainty in determination of 
the t x A17 value (target thickness times solid angle, as 
described above), 20% uncertainty associated with the- 
oretical analysis (uncertainty of the square of the ANC 
due to the variation of the optical potential parameters), 
10% uncertainty in the total resonance width and 10% 
uncertainty due to the initial a-d ANC. Thus, the S(0) 
value of l/2 + resonance determined in this experiment 



is (2.5 ± 0.7) x 10 6 MeVxb. The S(0) of the 6.356 MeV 
l/2 + state in 17 determined in this experiment is ten 
times smaller than that adopted in the NACRE compi- 
lation |3j and a factor of 5 larger than in [5j. The two 
channel R-matrix approach was used to calculate the to- 
tal S-factor of the 13 C(cv,n) process. All resonances from 
4.14 MeV (neutron decay threshold) to 8.2 MeV exci- 
tation energy in 17 were included. Parameters of the 
resonances were taken from [TT| . The a-reduced widths 
of resonances were fitted to reproduce the 13 C(a,n) ex- 
perimental data Q. It was found that it is necessary 
to introduce a constant 0.4 x 10 6 MeVxb non-resonance 
contribution to fit the experimental S factor data at the 
lowest energy region. The solid curve in Fig. |2| corre- 
sponds to the best R-matrix fit, obtained as explained 
above, with contribution from l/2 + state. Our calcu- 
lated reaction rates are shown in Figure 0] The best fit 
rate is shown as a solid curve in Fig. 0] in comparison 
with the NACRE adopted rate (dash-dotted curve). At 
temperatures above 0.3 GK, where contribution of l/2 + 
is small, it is identical to the curve adopted in NACRE, 
however, at temperatures which are most significant for 
the s-process in AGB stars, 0.08-0.1 GK, the reaction 
rate is smaller by a factor of 3 than that adopted in the 
NACRE compilation. Uncertainty in this astrophysically 
important reaction rate is now reduced from ^300% to 
15%. The inset in Fig. H shows the NACRE compila- 
tion adopted rate (dash-dotted curve) with uncertainty 
band (yellow) and the reaction rate obtained in this work 
(solid line) with uncertainty band shown in red. Numer- 
ical values of reaction rate for 0.08 - 0.1 GK temperature 
range are given in Table 1TT1 
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Table I: Parameters of Optical potentials used in DWBA calculations. 

Vo av r v W Ws a w r w r c V so a so r ao 

Channel Potential (MeV) (fm) (fm) (MeV) (MeV) (fm) (fm) (fm) (MeV) (fm) (fm) Ref. 

6 Li + 13 C LCI 134.0 0.68 1.50 - 11.1 0.68 1.50 1.50 - - [17] 

d + 17 DOl 105.0 0.86 1.02 - 15.0 0.65 1.42 1.40 6.0 0.86 1.02 [15J 

d + 13 C DC1 79.5 0.80 1.25 10.0 - 0.80 1.25 1.25 6.0 0.80 1.25 [18] 

i 

R = roA^ \ tq = 1.25 fm and a — 0.68 fm were used for a + d and a + 13 C form factor potentials with V fitted to reproduce 

binding energy. 



Table II: The rate of 13 C(a, n) reaction at temperatures from 
0.08 to 0.1 GK. The rate obtained in this work is compared 
with the rate published in the NACRE 3] compilation. Units 
are [cm 3 mol _1 s _1 ], exp stands for 10 exp . High and low values 
were calculated assuming 26% uncertainty of l/2 + 6.356 MeV 
resonance contribution. 
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reaction rate at stellar temperatures. Combination of 
the sub-Coulomb a-transfer reaction and application of 
the ANC technique in the analysis of experimental data 
practically eliminates all dependence of the results on 
model parameters, making this approach a very valuable 
tool for future studies of astrophysically important reac- 
tion rates with both stable and radioactive beams. The 
13 C(ev,n) reaction rate at stellar temperatures was found 
to be lower by a factor of 3 than previously adopted , 
also uncertainty in this reaction rate was greatly reduced. 
It would be of great interest to incorporate the new re- 
action rate into the s-process calculations in AGB stars. 



In summary, in this Letter we developed an indirect 
technique which allows measurement of the astrophysical 
S(0) factor of sub-threshold, particle unbound resonances 
and applied this technique to measure the contribution 
of the 1/2+ 6.356 MeV resonance in 17 to the 13 C(a,n) 
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